Long-distance entanglement in spin systems 
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Most quantum system with short-ranged interactions show a fast decay of entanglement with the 
distance. In this Letter, we focus on the peculiarity of some systems to distribute entanglement 
between distant parties. Even in realistic models, like the spin-1 Heisenberg chain, sizable entan- 
glement is present between arbitrarily distant particles. We show that long-distance entanglement 
appears for values of the microscopic parameters which do not coincide with known quantum critical 
points, hence signaling a transition detected only by genuine quantum correlations. 

PACS numbers: 03.67.Mn, 75.10.Pq, 05.70.Jk 



Entanglement generation and distribution is a prob- 
lem of central importance in performing quantum infor- 
mation (QI) tasks, like teleportation f]| and quantum 
cryptography [2|- Typically, the entanglement between 
parties is created by means of a direct interaction. Since 
entanglement needs the presence of strong correlations, 
low dimensional systems, as for example antiferromag- 
netic spin chains, offer a natural source of entanglement. 
In most systems with short-range interactions, the entan- 
glement between a pair of particles decays rapidly with 
the distance (generally even more rapidly than standard 
correlationsj. For example, in the Ising model with trans- 
verse field [3, |j| the concurrence vanishes for distances 
larger than two sites, while in the Heisenberg model (a,la| 
it is restricted only to nearest neighbors. 

Prom the QI perspective, it would be attracting to cre- 
ate sizable entanglement between particles that are lo- 
cated at a distance larger than a few sites. Along this di- 
rection, the localizable entanglement (LE) was conceived 
with the idea of exploiting spin chains as quantum chan- 
nels P, y|. The LE measures the average entanglement 
localized between a couple of distant points, after per- 
forming optimal local measurements onto the rest. 

In the present Letter, we show that already the ground 
state (GS) of various models widely used in condensed 
matter physics offer the possibility to entangle parties 
that are arbitrarily far apart. This fact, naturally leads 
to the concept of long distance entanglement (LDE) as 
a sort of quantum order parameter. As discussed in the 
following, the onset of LDE does not coincide with known 
quantum phase transitions (QPT's) of the systems we 
have examined. 

Let us consider two sites A and B that interact with a 
many-body system C. The distance d between A and B 
is set by the individual short-ranged interactions in the 
subsystem C (see figure below). 
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According to our definition, given a bipartite measure 
of entanglement E(p), we have LDE if 
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where pab — Trc |*) (^| is the reduced density matrix 
of the subsystem A and B, and |^) is the total wavefunc- 
tion. The introduction of two special points, or probes, 
is essential here since the property of monogamy J2| lim- 
its to two the number of particles maximally entangled. 
The basic idea comes from the observation that if we 
wish to locate a great amount of entanglement between 
two selected qubits, we are forced to exclude entangle- 
ment with the rest. Specifically, we have considered cases 
where A and B represent end spins in an open chain or 
additional spins (probes) that interact with selected sites 
in the chain. In condensed matter systems, these might 
be impurities, defects or even scattering particles |inl |. 

As a first criterion, we expect to have a nonvanish- 
ing LDE between A and B when their interactions with 
C are small compared to the typical interactions con- 
tained in C. Otherwise, A or B would develop too strong 
correlations with the closest degrees of freedom in C, ex- 
cluding the possibility to form LDE. On the other side, 
strong correlations among the particles in C will tend 
to avoid entanglement between C and the probes. In 
this sense, strongly correlated quantum systems, like an- 
tiferromagnetic spin systems, are good candidates to do 
the job. In particular, spin-1/2 antiferromagnetic sys- 
tems admit a simple picture based on resonating valence 
bonds (RVB) |ll|. If a state is a total spin singlet, then it 
may be approximated by all the possible RVB configura- 
tions, each one with a given weight. Resonances between 
various configurations destroy entanglement. The vari- 
ational idea for favoring a singlet between two selected 
sites (A and B) is to induce a large weight for all the 
RVB configurations that link pair of particles inside C 
by increasing here the interactions. 

In the following, we will present some mechanisms able 
to produce LDE in spin-1/2 and spin-1 chains. 



The dimerized-frustrated model. In the S = 1/2 an- 
tiferromagnetic isotropic Heisenberg chain, each spin is 
highly entangled with its nearest-neighbors [l3|- Instead 
we consider here the dimerized chain with frustration, 
well-known for its connections with spin-Peierls [l3| and 
ladder compounds, whose Hamiltonian is 
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where <j'^ , v = x,y,z are Pauli matrices. For 6 — 
the system is gapless up to ac ~ 0.241, where the GS 
spontaneously dimerizes and becomes doubly degenerate. 
In the Majumdar-Ghosh line 5 -\- 2a = 1 the system is 
made only by short-ranged singlets. 

We choose L even and open boundary conditions 
(OBC), with the aim to study the entanglement between 
the two spin-1/2 at the end points. First, let us look 
at two limit cases. For 5 = —1 and a — the GS is 
dimerized onto the "odd" bonds, 
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Figure 1: The finite size study on the concurrence shows the 
presence of LDE in model Q for 5 > Sria). Data were 
obtained keeping 256 DMRG states, with a truncation error 
smaller than 10"^". 



where the entanglement is localized in pairs of nearest- 
neighbors. More interesting for us is the case 5 = 1 and 
a = 0, where two spins are "left alone" as in the following 
figure 



and the GS is fourfold degenerate. The basic idea for 
concentrating a large amount of entanglement between 
end spins (A and B) is related to their tendency to 
form a global singlet in the GS for any 5^1. Hence 
the two end states are forced to develop strong corre- 
lations towards the formation of a long-distance singlet 
state I*") = (Iti) - lit)) /V^- This phenomenon can be 
thought of as a long range antiferromagnetic interaction 
mediated by the other spins in the chain. Accordingly, 
the states in Stot = 1 form a triplet of excitations. 

As a measure of entanglement, we adopt the concur- 
rence 14i. Given the SU(2) invariance of the GS, the 
spin-spin correlations ^^-^ = {a'^a'^) /4 are the same 
for every u. In addition, the magnetization is zero, so 
that the concurrence between A and B reduces simply to 
Cab = 2 max {0, 2 \Yab I ~ 7ab ~ 1/4}- The concurrence 
is nonzero if the antiferromagnetic correlations between 
A and B are sufficiently strong: 7^3 < —1/12. 

First, we have performed some numerical evaluations 
on the GS using the density matrix renormalization 
group (DMRG) method |l3|. The end-to-end concur- 
rence Cab is plotted in Fig. Q]as a function of the system 
size L for several values of 5 and a. The numerical data 
put in evidence the presence of LDE as well as the rapid 
achievement of the asymptotic value. This latter fea- 
ture is consistent with the small correlation length in the 
regime 5 > 0.10 (see e.g. jl3|) and allows us to study 



shorter chains by means of an exact diagonalization pro- 
gram based on the Lanczos method. 

The numerical results summarized in Fig.[2lshows that 
the end-to-end concurrence grows rapidly with 5 starting 
from a threshold value 5t{o)- For J = no LDE is 
generated, and this is related to the tendency of the first 
spin to entangle with the second, as found in Ref. |l(i |. 
This is also consistent with the absence of surface order 
in an open S" = 1/2 Heisenberg chain jl7|. 

The inclusion of a < tends to favor a classical Neel 
state, so it is expected to destroy entanglement. On the 
contrary, a > is seen to enhance the end-qubit concur- 
rence, as frustration favors quantum fluctuations. This is 
shown in the inset of Fig. [21 where (St (a) decreases with 
a, reaching a minimum for a ~ 0.5. In the limit |a| 3> 1 
the entanglement gets suppressed as the probes belong 
to two separated chains. Remarkably, from Figs.Q]and[2l 
it emerges that the entanglement grows with the system 
size L. 

Spin-1 Chain. An important class of spin-1 models 
is given by the Heisenberg chain with biquadratic inter- 
actions, 
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that has attracted much interest both for the study of 
hidden order |l^ and for optical lattice implementations 
0. At the AKLT point 3= 1/3, the GS is given by a 
valence bond solid (VBS) [l^, where each spin 1 is repre- 
sented by a couple of spin 1/2, provided the antisymmet- 
ric state is projected out. The VBS state is constructed 
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Figure 3: Entangled (white) and separable (shaded) states 
of two globally SU(2)-invariant qutrits are completely deter- 
mined by means of (SJSg) and {{SX) (S^) ). 



state to be entangled is that of having positive negativity 
12(11 ■ defined as 



Figure 2: Concurrence for the two end-qubit state as a func- 
tion of dimerization 5 for some values of a. Exact calculation 
on a chain of length L = 24. The LDE increases steeply above 
a threshold and is enhanced by frustration. Inset: threshold 
value of dimerization St (a), above which end-to-end concur- 
rence starts to be nonzero for lengths L — 12, 16, 20. The 
crosses are the infinite size extrapolations of DMRG data with 
L up to 100. 



by forming short-ranged singlets between nearest neigh- 
bor 5=1/2 states and then symmetrizing local pairs 
to get back 5* = 1 states. For OBC there remains free 
effective spin-1/2 particles at the endpoints responsible 
for a fourfold degeneracy, in an analogous way as in the 
the model (QJ with 6=1 and a = 0. Away from /3 = 1/3 
the degeneracy is lifted, the GS is a total singlet Stot = 
and other valence bond configurations give contribution 
to the GS. Anyway, the VBS state is still a good ap- 
proximation for a wide range of /3's, in particular at the 
Heisenberg point /3 = 0. Due to strong correlations in 
the bulk, the two 5* = 1/2 end spins tend to organize as 
a |^~) Bell state, in order to give rise to a total singlet. 

Two different measures were considered to quantify 
the entanglement between two spin-1. The VBS pic- 
ture suggests the definition of the partial concurrence 
(VC) as the amount of entanglement between the spin- 
1/2 belonging to different spin-1 particles. One advan- 
tage of the VC is that it depends only on the z-z cor- 
relator: VCa,b = 2max{0,2|77^^B|-r7X^B -1/4}, with 
ry^B = (S'^S'b)/4. The symmetrization procedure dis- 
tributes the entanglement among the four qubit, so that 
the maximal possible value of the VC is 1/2. However, 
the VC may fail in detecting genuine qutrit entangle- 
ment, which is generally hard to quantify. In fact, for 
qutrit mixed states there is no simple expression for the 
entanglement of formation nor a simple criterion for sep- 
arability is known. Nevertheless, for SU(2)-rotationally 
invariant states, a necessary and sufficient condition for a 
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where p^g stands for the partial transpose with re- 
spect to subsystem A and \\G\\-^ = trVCG^. Cal- 
culating the negativity of a general SU(2)-invariant 
state, parametrized by the quantities {S]^S^) and 
{(Sj^) (S'b) ), we are able to recognize the separable 
states. All the possible states fall inside the triangle in 
Fig. whereas the shaded area represents the separable 
states. 

In the AKLT case f3 = 1/3, we choose the singlet 
among the four degenerate states, because this is the 
state one would approach by letting /3 ^ 1/3 and cor- 
responds to the GS of the periodic chain. From the 
exact solution, one finds (SfSl) = —{{Sf) (S'£) ) ~ 



-4/9 l + 6(-l)'^e 



_nip-L/o 



where ^aklt = l/ln(3) 



is the bulk AKLT correlation length. It follows that in 
the thermodynamic limit VC = 1/6 and A/" = 2/9, where 
both values are approached exponentially fast. This con- 
firms the hypothesis that we have qubit as well as qutrit 
entanglement. 

At the Heisenberg point /? = with OBC it is well es- 
tablished the presence of surface order lhnL^oo{SfSl) = 
—0.28306484(1) [21| approached also in this case expo- 
nentially fast with a bulk correlation length ^h "^ Q- 
With accurate DMRG simulation up to 100 sites we 
could establish a similar behavior for the correlations 
{(Si) (SI) ) with asymptotic value very close to 4/9. 
These data imply the existence of LDE in the Heisenberg 
model detected by a nonzero negativity Af = 0.0608426, 
even if qubit entanglement vanishes, i.e. VC = 0. We note 
here that both the Heisenberg (H) and the AKLT points 

(see Fig. 13 lie on the line where {{Slf [Slf) ^^^ 
ii^i) ){{^l) ) — 4/9, which means that the nonzero 
spins (effective charges) are uncorrelated. Further en- 
hancement of the LDE may be achieved in spin-1 models 
that present also end-to-end charge correlations. 



S — 1/2 Heisenberg model with probes. So far, we 
have considered situations where the probes are located 
at the end points of a chain. Now we consider a different 
case: a Heisenberg chain of length L and two additional 
S = 1/2 probes ta and tq 

L 

i=i 

where r is a vector of Pauli matrices. The spin-probe A 
interacts with the site 1, while B is connected to the site 
d + 1. The correlations between A and B will depend 
only by their distance d + 2, having assumed periodic 
boundary conditions (PBC) with ct^,!.]^ = u". 

By means of exact diagonaHzation, we have studied 
the concurrence CAB(d), varying the distance at fixed L. 
The results are illustrated in Fig. ^ for a closed chain of 
L — 2Q (plus 2 probes). We have rejected even values 
of d, as in these cases the GS is threefold degenerate, 
belonging to the sector Stot — 1- Moreover, due to the 
PBC's the maximum distance is reached at half chain, 
d — 13. When Jp — 1, no probe entanglement is found, 
since CA^{d) = for every d. As we expected according 
to our considerations above. Cab is enhanced by weak- 
ening the interactions between the probes and the spin 
chain. Already for Jp = 0.3, the entanglement is nonzero 
for every (odd) value of d and remarkably at Jp = 0.1 the 
probes are almost completely entangled. Finite size scal- 
ing of the concurrence between maximally distant probes, 
C^aid — L/2, L) exhibits a slow decrease of the concur- 
rence with L and it remains an open question whether it 
survives at the thermodynamic limit. 

In addition, we just mention that changing the sign 
of the probe interactions to ferromagnetic, Jp < 0, the 
concurrence increases further, extending the possibility 
of tailoring interactions that yield efficient entanglement 
creation at large distance. A similar behavior is observed 
by placing the probes at the ends of an open chain. As 
above, the finite-size effects are non-negligible due to the 
critical nature of the bulk. Preliminary DMRG calcula- 
tions with L up to 100 leave open the possibility of having 
LDE for ferromagnetic probe interactions. 

Conclusions. With this Letter, we aim to bring to 
the attention of the QI community a large class of spin- 
1/2 and spin-1 models capable of creating entanglement 
between distant parties. On the one hand, this property 
opens up the possibility to engineer QI devices like en- 
tanglers and quantum channels using strongly correlated 
low-dimensional systems. In particular, the phenomenon 
of concentrating the entanglement on the border of finite- 
size system seems to be particularly suited for optical lat- 
tice simulations. On the other hand, we observe that the 
transition point where genuinely quantum correlations, 
signaled by the concurrence, extend to long distance does 
not coincide with known QPT's. How this issue embod- 
ies in the statistical mechanics framework is a challeng- 




Figure 4: Concurrence between probes attached to a Heisen- 
berg chain of length L = 26 as a function of the distance. The 
calculation was done for various values of Jp. A dramatic in- 
crease of entanglement between distant probes appears as Jp 
is lowered. 



ing question. Conversely, local measures of entanglement 
show a singular behavior at QPT's that comes from the 
most relevant operator |22]. Specifically, in this work 
we have considered models with SU(2) symmetry, which 
is common in nature and help to make the calculations 
easier. Nonetheless, we verified that the results regarding 
LDE apply also to non SU(2)-symmetric cases. Further 
work is in progress in order to extend the investigation 
of LDE on other models, including electronic systems. 

Many thanks to H. Nishimori for providing the Lanc- 
zos code TITPACK. We are grateful to E. Ercolessi, F. 
Ortolani, M. Paris and S. Pasini for useful discussions 
and especially to G. Morandi for a careful reading of 
the manuscript. This work was supported by the TMR 
network EUCLID (No. HPRN-CT-2002-00325), and the 
COFIN projects 2002024522_001 and 2003029498_013. 



[1] D. Boschi, S. Branca, F. D. Martini, L. Hardy, and 

S. Popescu, Phys. Rev. Lett. 80, 1121 (1998). 
[2] N. Gisin, G.Ribordy, W. Tittel, and H.Zbinden, Rev. 

Mod. Phys. 74, 145 (2002). 
[3] T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 

032110 (2002). 
[4] A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature 

(London) 416, 608 (2002). 
[5] M. Takahashi, J. Phys. C 10, 1289 (1977), cond- 

mat/9708087. 
[6] B.-Q. Jin and V. Korepin, Phys. Rev. A 69, 062314 

(2004). 
[7] F. Verstraete, M. Popp, and J. I. Cirac, Phys. Rev. Lett. 

92, 027901 (2004). 
[8] L. C. Venuti and M. Roncaglia, Phys. Rev. Lett. 94, 

207207 (2005). 
[9] V. Coffman, J. Kundu, and W. Wootters, Phys. Rev. A 



61, 052306 (2000). 

[10] G. D. Chiara, C. Brukner, R. Fazio, G. M. Palma, and 
V. Vedral, quant-ph/0505107. 

[11] P. W. Anderson, Science 235, 1196 (1987). 

[12] K. M. O'Connor and W. K. Wootters, Phys. Rev. A 63, 
052302 (2001). 

[13] D. Augier, D. Poilblanc, S. Haas, A. Delia, and 
E. Dagotto, Phys. Rev. B 56, R5732 (1997). 

[14] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[15] U. SchoUwock, Rev. Mod. Phys. 77, 259 (2005). 

[16] N. Laflorencie, E. S. Sorensen, M.-S. Chang, and I. Af- 
fleck, Phys. Rev. Lett. 96, 100603 (2006). 



[17] F. C. Alcaraz, M. Baake, U. Grimm, and V. Rittenberg, 

J. Phys. A 21, L117 (1988). 
[18] L Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, 

Comm. Math. Phys. 115, 477 (1988). 
[19] J. J. Garcia-RipoU, M. A. Martm-Delgado, and J. L 

Cirac, Phys. Rev. Lett. 93, 250405 (2004). 
[20] H.-P. Breuer, J. Phys. A: Math. Gen. 38, 9019 (2005). 
[21] G. Path, O. Legeza, P. Lajko, and F. Igloi, 

cond-mat/0603078. 
[22] L. C. Venuti, C. Degli Esposti Boschi, M. Roncaglia, and 

A. Scaramucci, Phys. Rev. A 73, 010303(R) (2006). 



